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We study the superfluid-Mott insulator (SF-MI) transition in an one-dimensional optical lattice 
system, and employ the Bose-Hubbard model in two dimension with a combined potential of an 
optical lattice in one direction and a confining harmonic trap in the other direction, which we 
refer to as the pseudo one-dimension Bose-Hubbard model. There some excited states with respect 
to the harmonic trap are considered. The Mott lobes shrink in the fj, and J directions of the 
[i-J phase diagram. The shrinkage occurs because the interactions involving the excited states 
become weaker than that between particles in the ground state. The dispersion of the in-site 
particle increases because the energy spacing between the eigenstates of the Hamiltonian decreases 
at finite temperature. The excited states significantly affect the robustness of the MI phase at finite 
temperate. 

PACS numbers: 67.85.Hj, 67.85.Bc, 03.75.Lm 
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I. INTRODUCTION 

Since the realization of Bose-Einstein condensates in 
the cold atomic systems |1| , many studies both from the- 
oretical and experimental sides have been done. Vari- 
ous parameters, say the strength of the inter-atomic in- 
teraction and the configuration of the confinement po- 
tential, are experimentally well-controllable in the cold 
atomic systems Q , which offer many intriguing phenom- 
ena such as the superfluid-Mott insulator (SF-MI) transi- 
tion. Moreover, it is possible to realize a low dimensional 
system simply by tightening the confinement potential in 
one or two directions, and it provides a good testing field 
to study a low dimensional physics. 

Experimentally, the SF-MI transition in this system, 
which is the quantum phase transition, has been observed 
when the depth of the lattice potential is increased Q, 
and it was followed by the observations in several systems 
such as the multi-component one [4, the multi-band one 
and the low-dimensional ones jg, [tJ- Theoretically, 
the strongly-correlated bosonic system in optical lattice 
is well described by the Bose-Hubbard model , accord- 
ing to which the possibility of the SF-MI transition was 
expected before its experimental observation [9\. This 
model has been analyzed in various situations, i.e., at 
finite temperature [ToL [Tlj , for the multi-component sys- 
tems [T^ - ll^ . for the multi-band ones [H, El* and for 
the low-dimensional one [l7| . 

The one-dimensional optical lattice system is real- 
ized experimentally in a combined potential of a one- 
dimensional optical lattice and a cigar-shaped harmonic 
trap. When the harmonic trap is tight in the directions 
perpendicular to the optical lattice and the optical lattice 
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itself is deep, this system can be described in a good ap- 
proximation by the one-dimensional (ID) Bose-Hubbard 
model, restricted to the first band. While the ID Bose- 
Hubbard model without higher bands has been investi- 
gated in some detail both at zero and finite temperatures 
11], its validity is not clear at finite temperature. It is 
reported that the contributions of the higher excitations 
of the harmonic trap cannot be neglected for the one- 
dimensional system with Bose-Einstein condensate at fi- 
nite temperature (l8| . 

Our aim of this paper is to investigate the temperature 
dependence of the SF-MI transition in a one-dimensional 
optical lattice system. To this end, we concentrate on a 
cold atomic system in two dimension with a combined po- 
tential of an optical lattice in x-direction and a confining 
harmonic trap in y-direction, and suppose that the opti- 
cal lattice is so deep that the ID Bose-Hubbard Hamilto- 
nian without higher bands describes the dynamics in re- 
direction, but that some excited states in y-direction have 
to be taken account of because the harmonic potential is 
not tight enough to forbid transitions into the excited 
states absolutely. We refer to this model as the pseudo 
one-dimensional (P1D) Bose-Hubbard model throughout 
this paper. The reason for analyzing the P1D model here 
is that while it is simple and easy to analyze, and it can 
be realized experimentally with a merit that the strength 
of the harmonic trap is variable. It is also pointed out 
that the structure of the transition terms in the P1D 
Bose-Hubbard Hamiltonian is similar to those in the spin 
transition [l3|, 0] or the one between the Bloch-bands 
[Tfil HH , implying that our model serves as a prototype 
to understand the physics of the low dimensional SF-MI 
transition in more general situations. 

We study the P1D Bose-Hubbard model both analyt- 
ically and numerically in this paper, and show how the 
transition terms affect the SF-MI phase diagram at zero 
and finite temperatures. The main result at zero tern- 
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perature is the shrinkage of the Mott lobes, which origi- 
nates from a decrease in the effective on-site interaction 
due to the transitions between the harmonic states in y- 
direction. One may expect that such deformations of the 
Mott lobes can be fitted simply by adjusting an effective 
on-site interaction in the ID model, but it is not true 
and they reflect the complex structure of interaction in 
the P1D model. Although the phase diagram for the P1D 
model reduces to that for the ID one in the tight confine- 
ment limit, we note that the shrinkage is non- negligible 
even for comparatively tight confinement. The contribu- 
tions of the transition in y-direction is more enhanced 
at finite temperature. The Mott lobes disappear gradu- 
ally with an increase in temperature due to the thermal 
fluctuation. While the disappearance in the ID model 
occurs independently of the value of the chemical poten- 
tial , it is not so in the P1D model. This is explained as 
follows: the effective on-site interaction becomes weaker 
for particles in the excited states, and the fraction of the 
excited particles becomes larger as the chemical poten- 
tial is increased, which implies that the hopping term 
becomes relatively significant for the larger chemical po- 
tential. This way we have the dependence of the chemical 
potential in the phase diagram, calculated from the P1D 
model. We have confirmed through our analysis that it 
is only at very low density and temperature that the ID 
Bose-Hubbard model without higher bands and the P1D 
one give almost the same results, and therefore that the 
validity of the ID model is restrictive. 

This paper is organized as follows. In Sec. (Til we intro- 
duce the P1D Bose-Hubbard model. The phase diagram 
in the zero hopping limit is obtained analytically. In 
Sec. Mil we show the phase diagram for non-zero hop- 
ping both at zero and finite temperatures by employing 
the mean-field approximation. Section IIVI is devoted to 
summary and discussion. 



II. PSEUDO ONE-DIMENSIONAL 
BOSE-HUBBARD MODEL 



model throughout this paper for simplicity, an extension 
to three-dimensional systems, e.g. one with a combined 
potential of an optical lattice in x-direction and a confin- 
ing harmonic trap in j/z-plane, is straightforward. 
The field operator ip(x,y) is expanded as 



(2) 



in the complete set of wave functions of harmonic oscil- 
lation, 



H 2 d 2 

1T--T2 + V ^(v) ) u t(u) = E e ue(y) , (3) 
Ira dy z ' 



where E t = fkj y (£ + 1/2) (I = 0, 1,2, ■••). Then we 
expand the field operator ipe( x ) m t ne complete set of 
the Wannier functions Wi(x) with site indices i as 



(4) 



We assume here and hereafter that the optical lattice 
potential is so deep that only the states in the lowest 
band contribute, but take account of the excited states 
(£ 0) of harmonic oscillation in y-direction. Under the 
tight-binding approximation we obtain the Hamiltonian, 
which we refer to as the P1D Bose-Hubbard Hamiltonian, 
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We start with the following Hamiltonian to describe 
the cold Bose atomic gas system in two dimension with 
a combined potential of an optical lattice in ir-direction 
and a confining harmonic trap in y-direction as follows: 

f r - / h 2 o 2 \ - 

H= Idxdy Tl)^x,y)[ - —— + V opt {x) jip{x,y) 



h 2 d 2 

+ft {x,y)[ -^q~2 + IWiO -**)''< ' ■ !' 1 
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ft (x, y)ft (x, y)ftx, y)ftx, y) 



(1) 



where V ap t{x) t Vhar(y) 



;y 2 /2, /i, and g represent 



the optical lattice potential, the harmonic potential, the 
chemical potential, and the coupling constant, respec- 
tively. Note that no harmonic potential in a;-direction is 
considered. Although we study only this two-dimensional 



where 



J 



f ( ti 2 d 2 \ 
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gtilttsli = 9 J dy u* £l (y)u; 2 (y)u ea (y)u ei (y) , 
UlililaU = 9I1I2I3U dxw*(x)w*(x)wi(x)wi(x) . 



(9) 
(10) 

(11) 

4 de- 
pends on £ (the quantum number of the harmonic oscil- 
lation) while the hopping coefficient J is independent of 
I. In Table [l] the ratios of the on-site interaction coef- 
ficients to the ground state interaction are shown. All 
the on-site interaction coefficients involving the excited 
states are smaller than the interaction coefficient C/qooo- 



Note that the on-site interaction coefficient Ut x i 
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TABLE I. The ratios of the on-site interaction coefficients 
Ui 1 e 2 e 3 t 4 to the ground state interaction IToooo-The coefficients 
have the symmetric property, for example, (Toon = 
C^oioi — C^noo, etc- 
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The harmonic frequency io y representing the strength 
of the harmonic potential comes into the dynamics of 
our model in two ways, namely through the energy-level 
spacing AE y — — Ei = fno y and the on-site inter- 
action L^ 1 £ 2 £ 3 £ 4 . The latter dependence is expressed as 
Ui 1 i 2 i 3 i A r-j y/LJy~ universally for any £, and is absorbed 
into rescaling the phase diagram by using the parameters 
/V^oooo an( i J/Uoooo- Note that it does not mean that 
the P1D model reduces to the ID model by such rescal- 
ing. We pay little attention to this trivial uj y dependence, 
but focus on the nontrivial effect of the transitions be- 
tween the harmonic states, which is absent in the ID 
model. 

We assume that the confinement harmonic potential 
is relatively tight but is not so tight as all the excited 
states can be neglected. Then, the infinite sum over £ 
can be approximated by a finite sum of some lower energy 
states. Because of the even-odd selection rule for the on- 



site interaction coefficient Vi 1 i 2 i 3 i i , the transition of one 
particle from the ground state to the first excited one is 
forbidden. The transition with the lowest energy from 
the ground state is either the simultaneous one of two 
particles to the first exited state or the one of one particle 
to the second excited state. Hence, in order to estimate 
effects of the excited states at a minimum, we have to 
sum up the states not with £ = 0, 1 but with £ = 0, 1, 2 
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at least, and will take this minimal sum throughout this 
paper. 

Before full numerical calculations, let us turn to the 
zero hopping limit, J = 0. Then the hopping Hamilto- 
nian if hop in Eq. (J5J) is absent, and the total Hamilto- 
nian H — on-site is entirely separable with respect to 
the site index i. Note that hi in -f/on-sitc commutates 
with the in-site particle number operator at the i-site, 

hi = J2e ®t i®t,i ' an d hence the ground state of the total 
system is obtained simply by the diagonalization with a 
fixed in-site particle number. Hereafter, the trivial sub- 
script i is omitted. 

Explicitly we set up the following eigenequation for 
each fixed in-site particle number n and seek the lowest 
energy eigenstate, 

h\^) n = E n \^) n , (12) 
where the state is expanded as 

|*)n = X! ^V»[i+"i+"2 9i n o, n l> n 2)\n , n i, n 2j , 
n «l "2 

(13) 

with the direct product of the particle number states \ni) 
with £ = 0,1,2. For illustration we write down the matrix 
eigenequations for n — 1 and n = 2: 
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(16) 
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with 



$ 2 



/ff(2,0,0)\ 
3(1,0,1) 
3(0,0,2) 
3(0,2,0) 
3(1,1,0) 

Wo, 1,1)/ 



(17) 



Equation (1X1)) implies that the lowest eigenvalue of Ei 
is equal to £a )m i n = — /i. One can calculate the the low- 
est eigenvalue of E 2 , denoted by £?2,min, from Eq. (fT6|) . 
Similarly, the lowest energy and its eigenstate for each 
n are obtained. The results allow us to draw the phase 
diagram as indicated in Fig. [T] There the boundary be- 
tween n = 1 and n = 2 regions, for example, is a line on 
which Ei, min = E 2 ,mm- 
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TABLE II. The squared coefficient \g(n ,n 1 ,n 2 )\ 2 for AE y /U ooo = 3.0. 



in-site particle number 


| 3 (n,0,0)| 2 


| 3 (n-l,0,l)| 2 


| 3 (n-2,2,0)| 2 


others 


n = 2 


9.84 x HT 1 


7.88 x 10" a 


7.47 x 1Q-' 6 


8.37 x 10" 4 


n = 3 


9.35 x lO" 1 


4.23 x 10~ 2 


2.11 x 10" 2 


1.49 x 10" 3 


n = 4 


8.57 x 10" 1 


1.07 x 10" 1 


3.48 x 10~ 2 


1.15 x 10~ 3 



Obviously, the dispersion of the in-site particle number 
is equal to zero as a result of its conservation, and the 
state of the whole system is in a Mott phase. But the 
particle number in each £-state can fluctuate. The dis- 
persions of the particle number in each ^-state are shown 
in Fig. [5J where they varies stepwise as functions of /i. 
For comparison, we recall the corresponding ID Bose- 
Hubbard model, 

#1D = + + X! ^' 1D ' ( 18 ) 



hi,iT> = -fJ-al&i + -C/qooo alalaiai . 



(19) 



In the limit of J = 0, we consider the eigenequation 
hii)\^)n — E n \^f) n for each fixed n, and its lowest eigen- 
value is found to have a very simple form of 



E,, 



n(n — 1) 
-n/x H C/qooo ■ 



(20) 



which implies that the phase boundaries are at /i/J/oooo = 
1, 2, • • • [Hj. We can see from Fig. Q]that the existence 
of the excited states £ = 1,2 in the P1D model pushes 
down the phase boundaries towards lower fi. It comes 
from the fact that any on-site interaction involving the 



excited states Ug 1 t 



is smaller than the interaction 



between ground state particles t/oooo- 

In Table UH we show the squared coefficients 
\g(no, n\, 722) | 2 for the in-cite particle numbers n = 2,3, 4. 
The probability of the ground state \n, 0, 0) is the greatest 
for each n, and the next leading probabilities, |n— 1, 0, 1) 
and \n— 2, 2, 0), are nonnegligible. On the other hand, the 
value in the column "others" , corresponding to a sum of 




FIG. 1. The phase diagram in the zero hopping limit, J = 0, 
for the P1D Bose-Hubbard model. The values of /i/f/oooo on 
the phase boundary line approach in the tight confinement 
limit, AEy/Uoooo — > 00 , to the ones for the ID Bose-Hubbard 
model, which are /z/t/oooo = 1, 2, • • • . 
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FIG. 2. (Color online) The dispersion of the particle number 
in the ground, first excited, and second excited states £ = 
0,1,2 for the P1D Bose-Hubbard model with AE y /Uoooo = 
3.0. 



contributions of the states with higher energies, is negli- 
gibly small within our parameter region. This guarantees 
the validity of our restriction to the energy levels up to 
£ = 2 when the energy-spacing AE^/C/oooo is relatively 
large. 



III. MEAN-FIELD APPROXIMATION AND 
NUMERICAL RESULTS 

In this section, we show numerical results of the SF-MI 
transition and the /i- J phase diagram for the P1D Bose- 
Hubbard model with the hopping Hamiltonian (the non- 
vanishing J) at zero and finite temperatures. To deal 
with the hopping Hamiltonian, we resort to the mean- 
field approximation [111 ] . 

Let us introduce the order parameter 



{at) 



Tr[e-^} 



(21) 



where (3 is the inverse temperature /3 = 1/ksT with 
the Boltzmann constant fee ■ Neglecting the fluctuation 
terms of the second order, we have an approximate ex- 
pression of 



",.,"'„/ - ®t,M,j + $l,jal,i ~ ®Xi®t,j ■ (22) 
Then, the approximate hopping Hamiltonian becomes 



(M) 
hop 



-j 



£[(* 



e,i+i 



i+l 



(23) 
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FIG. 3. (Color online) (a) The dispersion of the in-site par- 
ticle number for the ID model, (b) The dispersion of the 
in-site particle number, (c) That of the ground state particle 
number, (d) The absolute square of the ground state order 
parameter |$o,i| 2 for the P1D model with AE y /Uoooo = 3.0 
at zero temperature. The black colored regions where the dis- 
persions are zero in (a) and (b) denote the Mott lobes, and 
the n-th lobe from the bottom represents the Mi-phase with 
n in-site particles. The regions where the order parameter 
vanishes, as depicted in (d), correspond to the Mott lobes in 
(b). 

and the total Hamiltonian under the mean-field approxi- 
mation is given by ij( M ) = H^+H on -site ■ We calculate 
the lowest energy state by diagonalizing the total Hamil- 
tonian, when the order parameter $^ is determined self- 
consistently by minimizing the energy at zero tempera- 
ture or the free energy 

F = -/3~ 1 In Z(/3), (24) 

Z(/3)=Tr[e-^ <M) ], (25) 

at finite temperature. In this calculation, we consider the 
excited states up to the second excited state as stated in 

Sec. nu 



A. Results at Zero Temperature 

In Fig. [31 we show the phase diagrams for the ID model 
and for the P1D model with the energy-level spacing 
AEy/Uaoao = 3.0 at the zero temperature. Comparing 
the result for the P1D model with that for the ID model 
in Fig. |3J-(a), one can find the following two features. 
First, the Mott lobes for the P1D model shrink in the \x 
direction in Fig. |3j-(b). The rate of the shrinkage of the 
Mott lobes for the P1D model increase as the in-site par- 
ticle number increases. Secondly, the dispersion of the 
ground state particle number in the Mott lobes for the 
P1D model is nonzero and also increases as the in-site 
particle number increases, as in Fig. Er( c )- The order 
parameter does not necessarily become nonzero even in 
the nonzero dispersion region of the ground state particle 




FIG. 4. (Color online) (a) The dispersion of the in-site particle 
number and (b) the absolute square of the order parameter 
with /3(7oooo = 20.0 for the ID model. Figures (c) and (d) 
represent those with /3t/oooo = 8.0 , respectively. 

number in the Mott lobes, as in Figs. [3]-(c) and (d). The 
nonzero dispersion of the ground state particle number 
shows that the particle number of each state is not fixed 
and can fluctuate, but only the in-site particle number 
is fixed. This is attributed to the decrease of the ef- 
fective on-site interaction because the on-site interaction 
coefficients are smaller than the ground state interaction 
coefficient E/oooo- The numerical results show that the 
analysis in Sec. [TT| can be extended even to the case of 

We see from Figs. |3]-(a) and (b) that the Mott lobes 
for the P1D model also shrink in the J direction. The 
reason for the shrinkage in the J direction is as follows: 
the system is favorable to the SF-phase when the ratio of 
the hopping to the on-site interaction is large. The effec- 
tive interaction is smaller than that of the ground state 
particles, which means that the particles in the excited 
states can easily be in the SF-phase than the particles in 
the ground state because then the effect of the hopping 
is relatively large compared to the on-site interaction. 
As a result, in comparison with the system for the ID 
model, the systems for the P1D model is favorable to the 
SF-phase. 

We emphasize that the P1D model is essential in the 
above considerations since the particles in the excited 
states play a crucial role. 

B. Results at Finite Temperature 

In the previous subsection, we found that the effects 
of the excited states become more conspicuous for larger 
in-site particle numbers. It is therefore anticipated that 
the excited state is more important at finite tempera- 
ture. In this subsection, we analyze the effect of finite 
temperature for the P1D model. 

In Fig. 21 we show the particle number dispersions and 
the order parameters for the ID model at temperatures 
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FIG. 5. (Color online) (a) The dispersion of the in-site particle 
number and (b) the absolute square of the order parameter 
with l3U 000 o = 20.0 for the P1D model. Figures (c) and (d) 
represent those with /3J7oooo = 8.0 , respectively. 



P C^oooo = 20.0 and /3 E/oooo = 8.0. It is seen from Fig. H 
that the Mott lobes do not fit the region of the vanishing 
order parameter, the latter is wider than the former. We 
have some regions where the dispersion is non-vanishing 
but the order parameter vanishes, $o = 0, and interpret 
that they are in the normal-liquid phase 

The phase diagram for the P1D model is shown in 
Fig. [5j There the order parameter is that for the ground 
state. As in the ID model, the vanishing order parameter 
region covers the vanishing dispersion one. The normal- 
liquid region expands as temperature goes up. In a sim- 
ilar fashion to the zero temperature case, the shrinkage 
of the Mott lobes for the P1D model also is found at fi- 
nite temperature. The shape of the lobes for the order 
parameter in Fig. [5]-(b) and (d) is lost at finite temper- 
ature more clearly than for the ID model in Fig. HJ-(b) 
and (d). 

To investigate the temperature dependence of the dis- 
persion in more detail, we show the dispersions for small 
J, namely for J = in Fig. [5] and for J/Uqqoo = 0.02 
in Fig. [7] Comparing the results for the P1D model 
in Figs. 0(b) and GO- ( c ) with those for the ID model in 
Figs. |6]-(a) and [3(a), we see that the dispersion for the 
P1D model is larger than that for the ID model. In ad- 
dition, the dispersion for the ID model is independent 
of the chemical potential /Lt, but that for the P1D model 
depends on \i in such a manner that it is larger at larger 
/i. The order parameter for the P1D model begins to be 
non- vanishing in smaller fi region than for the ID model, 
as in Fig.[B]-(b) and (d). At higher temperature, the order 
parameter becomes smaller and the region in the normal- 
liquid phase expands equally for the ID and P1D models 
in Fig. El 

Figure [8J- (a) shows the dispersions of the average in- 
site particle number with (n) — 2,3,4 and J = 0. In 
both of the ID and P1D models, the dispersions increase 
as the temperature rises, but the increasing rate of the 
dispersion for the P1D model depends on the average 





FIG. 6. (Color online) The dispersions of the in-site particle 
numbers in the zero hopping limits J = for (a) the ID and 
(b) the P1D models with AE y /U 00 oo = 3.0 . 



00 





FIG. 7. (Color online) The dispersions of the in-site particle 
numbers (a) and the absolute square of the order parameter 
(b) with J/Uoooo — 0.02 for the ID model. Figures (c) and 
(d) represent those for the P1D model respectively. 



particle number (n). We also illustrate the temperature 
dependence of the increasing rates of the excited particle 
numbers for £ = 1, 2 in Fig.[5]-(b) . One can find in Fig.[S] 
that the increase in the excited particle numbers synchro- 
nizes with that in the dispersions. This suggests that the 
excited particles contribute to the large dispersions in the 
P1D model. 



TABLE III. In this table, E^ D) and E^ LLJ > represent the ra- 
th eigenvalues of the Hamiltonian for the ID and P1D systems 
with (n) = 3 and /3[/oooo — 10. These values are respectively 
the eigenvalues of the grand canonical Hamiltonian H — n/j,, so 
the values can be negative. The symbol A_E m _i jln represents 
the spacing between energy levels, A_E m _i im = E m — E m -\ . 
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-4.00 
0.500 
-3.19 
0.359 


-4.00 
0.00 

-3.19 
0.00 
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1.50 
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FIG. 8. (Color online) (a) The dispersions of the average 
in-site particle number with (h) = 2,3,4 and AE y /Uoooo = 
3.0 for the P1D and ID models. The dispersion for the ID 
model is independent of the average in-site particle number 
(n) . (b) The ratios of the excited particle numbers at finite 
temperature rii to those at zero temperature rii t T=o in the 
average in-site particle number (n) = 3 for the P1D model 
with AE y /U ooo = 3.0. 

To clarify the rapid increase of the dispersions for the 
P1D model, we investigate the eigenvalues of the Hamil- 
tonian, 

H^\V m ) =E m \* m ), (26) 

for the ID and P1D models. In Table. ITlTI the eigenval- 
ues of the Hamiltonian and the energy-level spacings are 
shown. One finds that the energy-level spacing for the 
P1D model AE^^ m is narrower than that for the ID 

model AE^\ m . The narrow energy- level spacing en- 
hances the contributions of the higher energy eigenstates 
m = 1,2, • • • at finite temperature, so that the disper- 
sions become large then. In addition, it turns out for 
(h) = 3 that the eigenstates belonging to the eigenvalues 
Ei and Ei are those with the in-site particle numbers 
ri = 4 and n = 2, respectively, while the eigenstate to 
belonging to Eq is that with n — 3 . The number of 
the excited particles grows as the in-site particle num- 
ber does according to Table [TTJ Thus the increase of the 
excited particles in Fig.[8]-(b) reflects the contribution of 
the eigenstate with Ei, equivalently of the state of the 
in-site particle number n = 4. 

The narrow energy-level spacing arises from the fact 
that the on-site interaction coefficients involving the ex- 
cited particles are small. This is explained in the follow- 
ing way. As the fraction of the excited particle number 
becomes large at the large average in-site particle num- 
ber (h) , the effective repulsive interaction becomes small 
and so are the eigenvalues E m (m = 1, 2, • • • ) . 

The above discussions indicate that the presence of the 
excited states affects the robustness of the Mott phase 
definitely and essentially. Although the approximation of 



the ID model is valid for very low density temperature, 
the effects of the excited states can not be neglected for 
describing more general situations. 

IV. SUMMARY AND DISCUSSION 

In this paper, we analyzed SF-MI transition at zero 
and finite temperature for the P1D Bose-Hubbard model. 
The P1D model, which differs from the ID model, con- 
tains the transition terms, which are similar to the 
spin transition and the transition between the Bloch 
bands. Specifically, the form of the multiband Bose- 
Hubbard Hamiltonian corresponds to that of the P1D 
Bose-Hubbard Hamiltonian, aside from the values of the 
on-site interaction coefficients U£ 1 i 2 i 3 i 4 and the hopping 
coefficients. We expect that our analysis is also useful in 
analyzing the multiband model. 

In zero temperature analysis, the SF-MI phase diagram 
for the P1D model is different from that for the ID model 
even for the relatively tight confinement. The major dif- 
ference is the shrinkage of the Mott lobes in the fj, direc- 
tion. The effect is attributed to decreasing the effective 
on-site interaction by the transition between the ground 
and excited states. The transition produces the fluctu- 
ation of the particle number in each £-state in the Mott 
lobes. The decrease of the on-site interaction enhances 
the hopping relatively, which leads to the shrinkage of the 
Mott lobes in the J direction as that in the /i direction. 
We note that the ID model can not take account of these 
effects by adjusting the parameters. 

At finite temperature, the dispersion of the in-site 
particle number increases as the temperature becomes 
higher. In addition, the region of the normal-liquid phase 
surrounding the Mott lobes expands. These effects occur 
both for the ID and P1D models due to the thermal 
fluctuation. The result for the ID model corresponds to 
Ref. [ll[. We note that the increase of the dispersion 
for the P1D model is enhanced more than that for the 
ID model. The tendency is more noticeable at higher 
temperature due to the increase of thermal fluctuations. 
These results are attributed to the weaker on-site inter- 
action coefficients involving the excited states. 

It have been seen that the robustness of the MI phase 
is affected at finite temperature by the presence of the 
excited particles. From the view of the local density 
approximation in Fig. [5], the P1D Bose-Hubbard model 
predicts that the Mott shell collapses more quickly than 
the ID model does. The validity of the ID model is re- 
stricted to very low density and temperature for realistic 
cold atomic gases. At high density and/or temperature, 
one has to to consider the excited states coming from the 
confinement potential. 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. [3] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

Wicman, and E. A. Cornell, Science 269, 198 (1995). I. Bloch, Nature (London) 415, 39 (2002). 

[2] I. Bloch, Nature Phys. 1, 23 (2005). [4] J. Catani, L. De Sarlo, G. Barontini, F. Minardi, and M. 



8 



Inguscio, Phys. Rev. A 77, 011603 (2008). 
[5] T. Miiller, S. Foiling, A. Widera, and I. Bloch, Phys. 

Rev. Lett. 99, 200405 (2007). 
[6] I. B. Spielman, W. D. Phillips, and J. V. Porto, Phys. 

Rev. Lett. 98, 080404 (2007). 
[7] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. 

Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 
[8] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. 

S. Fisher, Phys. Rev. B 40, 546 (1989). 
[9] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[10] K. Sheshadri, H. R. Krishnamurthy, R. Pandit, and T. 

V. Ramakrishnan, Europhys. Lett. 22, 257 (1993). 
[11] X. Lu and Y. Yu, Phys. Rev. A 74, 063615 (2006). 
[12] P. Buonsante, S. M. Giampaolo, F. Illuminati, V. Penna, 

and A. Vezzani, Phys. Rev. Lett. 100, 240402 (2008). 



[13] S. Tsuchiya, S. Kurihara, and Takashi Kimura, Phys. 

Rev. A 70, 043628 (2004). 
[14] M. Yamashita and M. W. Jack, Phys. Rev. A 76, 023606 

(2007) . 

[15] J. Larson, A. Collin, and J. -P. Martikainen, Phys. Rev. 

A 79 033603 (2009). 
[16] A. Collin, J. Larson, and J. -P. Martikainen, Phys. Rev. 

A 81, 023605 (2010). 
[17] N. Elstner and H. Monien, Phys. Rev. B 59, 12184 

(1999). 

[18] E. Arahata and T. Nikuni, Phys. Rev. A 77, 033610 

(2008) . 

[19] C. J. Pethick and H. Smith, Bose-Einstein Condensation 
in Dilute Gases Second Edition, (Cambridge University 
Press, 2008). 



